Network-like propagation of cell-level stress in random foams 
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Quasistatic simple shearing flow of random monodisperse soap froth is investigated by analyz- 
ing Surface Evolver simulations of spatially periodic foams. Elastic-plastic behavior is caused by 
irreversible topological rearrangements (Tls) that occur when Plateau's laws are violated; the first 
Tls occur at the elastic limit and at large strains frequent cascades of Tls, composed of one or 
more individual Tls, sustain the yield-stress plateau. The stress and shape anisotropy of individual 
cells is quantified by Q, a scalar measure derived from the interface tensor that gauges each cell's 
contribution to the global stress. During each Tl cascade, the connected set of cells with decreas- 
ing Q, called the stress release domain^ is network-like and highly non-local. Geometrically, the 
network-like nature of the stress release domains is corroborated through morphological analysis 
using the Euler characteristic. The stress release domain is distinctly different from the set of cells 
that change topology during a Tl cascade. Our results highlight the unique rheological behavior 
of foams, where complex large-scale cooperative rearrangements of foam cells are observed as a 
consequence of distinctly local events. 



Foams are complex fluids with a shear modulus and 
a yield stress that are linked to shear-induced deforma- 
tions of jammed nonspherical bubbles [U |2]. Shaving 
foam furnishes a common but illustrative rheological ex- 
perience; it is a soft solid at low stress yet flows at high 
stress. The defining mechanisms behind yield stress flu- 
ids are predominant under quasistatic conditions where 
the deformation rate is small and viscous forces are neg- 
ligible. The shear modulus and yield stress characterize 
the elastic-plastic response, and both increase as the vol- 
ume fraction of gas increases, reaching a maximum in 
the dry foam limit where the liquid volume fraction ap- 
proaches zero [H H]. Remarkably, foam is stiffest and 
strongest when it contains mostly gas. 

Gonsider a dry soap froth subjected to quasistatic 
shear (Figure [T]). The idealized microstructure consists 
of a network of surfaces that divide space into polyhedral 
cells and represent thin liquid films with negligible thick- 
ness. Dry foams in equilibrium satisfy Plateaus laws, 
which state that three films meet along each edge at 120°, 
and four edges meet at each vertex [5j. Plateau's laws 
govern the local geometry of the equilibrium film net- 
work as a consequence of surface energy (area) minimiza- 
tion, and therefore, determine the rheology of the soap 
froth. The response of foam to deformation is elastic until 
Plateau's laws are violated; typically, a cell edge length 
goes to zero causing more than four edges to meet at the 
merged cell vertex. This triggers an abrupt Tl topologi- 
cal transition (Figure |2|. The shape changes that result 
from the Tl can cause other edge lengths to vanish, and 
trigger additional Tls; the Tl cascade terminates when 
Plateau's Laws are fulfilled and equilibrium is restored. 
The instantaneous rearrangement causes a drop disconti- 
nuity in the macroscopic stress. Subsequent deformation 
is elastic until the occurrence of further Tl cascades, etc. 

Elastic response punctuated by discrete cell-neighbor 
switching is a hallmark of the simple 2D Princen model 




FIG. 1. (Golor online) A dry liquid foam with 512 cells per 
periodic unit cell at strains of and 0.6 respectively. The cells 
are colored by Q, a scalar measure of cell distortion, where 
blue through red symbolizes values through 1. The average 
Q over all cells are 0.244 and 0.425 respectively. 




FIG. 2. (Golor online) Illustration of a Tl topological tran- 
sition in three dimensions, during which cells 1 and 2 lose 
contact and the other two cells become neighbors. 



of a liquid honeycomb in simple shear [6 ; and also has 
been analyzed in 3D for the Kelvin foam [7] and Weaire- 
Phelan foam [8l under shear [9l [10]. Small systems, 
however, are incapable of exhibiting a stress plateau in 
2D [11 or in 3D, because of constraints due to spatial pe- 
riodicity and flow kinematics. Simulations of disordered 
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2D foams [2j [121 US] do exhibit a stress plateau, which 
is a piecewise continuous function that fluctuates about 
its average value. The stress for disordered 2D foams 
increases during small elastic deformations and is par- 
tially relieved by drop discontinuities caused by Tl cas- 
cades. Large cascades can result in strain localization, 
also called shear banding [14 . 

Experimental studies have addressed shear banding in 
quasi-2D foams, which consist of bubbles confined to a 
single layer in various ways p!5HT8] . The shear bands ap- 
pear to be caused by stress inhomogeneities or by viscous 
drag from confining glass plates [19]. Shear bands have 
also been found in simulations of 2D foams in straight 
channels by analyzing spatial and temporal correlations 
of topological transitions j2QH23] : however, simulations of 
fully periodic (unbounded) 2D foams under homogeneous 
shear have not been reported. 

Diffusing wave spectroscopy has been used to detect in- 
termittent shear-induced bubble rearrangements during 
steady flow of shaving foam [24]426] . Topological tran- 
sitions have also been directly observed [27] in a soap 
froth during steady shearing flow. Experiments have yet 
to detect shear banding in dry 3D foams |28H3Q] . 

A method for simulating the equilibrium microstruc- 
ture of random soap froth with controlled cell-size distri- 
butions j3]TF33^ employs molecular dynamics to generate 
dense packings of rigid spheres, which are used to produce 
Laguerre (weighted- Voronoi) tessellations. The tessella- 
tions are initial conditions for the Surface Evolver [34^, 
which minimizes surface area given cell volume con- 
straints. Convergence to a local surface area minimum 
requires numerous Tls that are triggered by cell edges 
going to zero length. Annealing via large-deformation, 
tension-compression cycles further reduces surface area 
and the number of short cell edges. Finally, if isotropic 
stress is desired, a slight distortion (recoil) of the cu- 
bic unit cell is required. The topological statistics of 
the resulting random monodisperse foams are in excel- 
lent agreement with Matzke's seminal experiments [35] , 

Quasistatic simple shearing flow is implemented in the 
Surface Evolver by repeatedly subjecting a periodic foam 
to small strain steps (e.g., A7 = 0.005), which involve 
affine deformation followed by equilibration to a state 
that once again satisfies Plateau's laws. If any cell edge 
length is less than some cutoff, that edge is eliminated 
and topological transitions are implemented until equi- 
librium is achieved. The cutoff is typically between 1% 
and 5% of the average edge length in the undeformed 
structure (where the total shear strain 7 = 0). 

The non-isotropic part of the effective macroscopic 
stress Tij of a dry soap froth [31 is evaluated as 
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where S is the set of all surfaces that divide the foam into 
cells. Here, a is the surface tension; Vf is the volume of 



FIG. 3. The foam energy E (top), shear stress r (middle), and 
normal stress differences Ni and A^2 (bottom) , normalized by 
surface tension and volume aV~^, plotted as a function of 
the total shear strain 7. Each plot contains two simulations 
with 216 (6^) cells and one with 512 (8^) cells. For these 
simulations, the initial unit cell is cubic, and therefore, r, Ni 
and N2 are non-zero at 7 = 0. 



the foam; 5ij is the Kronecker delta; is a local unit 
vector normal to the surface S; and ds is the differential 
area element. The factor of 2 occurs because each surface 
(soap film) has two sides. 

The foam energy density is defined as E = cjSf/Vf^ 
where Sf is the total surface area of the foam (counting 
both sides of each film). Rheological quantities of interest 
include the shear stress (r = T12) and the first and second 
normal stress differences {Ni = rn — r22, N2 = T22 — T33), 
where flow in the 1-direction depends on the 2-direction. 
Figure [3] shows the dependence of the energy and other 
rheological functions on the total strain 7. 

The foam response is initially elastic with r ~ 7 and 
£^,7Vi,A^2 ~ 7^- Eventually, all quantities are expected 
to reach a plateau (to within fluctuations). The strain 
at which the first Tl occurs is system dependent and de- 
fines the elastic limit and the onset of irreversibility. The 
elastic limit occurs at smaller strains and lower stress 
than the plateau, and there is some evidence of a shear 
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FIG. 4. (Left) A typical stress release domain for a purely elastic strain step consists of a few (normally between 3 and 10) 
disconnected cells. (Center) A typical plastic zone is composed of all cells that change topology across a topological transition. 
(Right) A typical stress release domain across a topological transition can contain hundreds of cells, and spans the entire foam. 
The foam has 512 cells and the strain is approximately 0.5. 



stress overshoot, which would correspond to a static yield 
stress. In the stress plateau regime, the stress is a piece- 
wise continuous function, composed of small elastic defor- 
mations followed by drop discontinuities caused by Tls. 

The normal stress differences reach a plateau more 
gradually than the shear stress, suggesting greater sen- 
sitivity to the foam microstructure. The ratio of — ^^"2 
to Ni in the initial elastic regime is slightly less than 1, 
comparable to the | predicted in [36^ . This ratio is pecu- 
liar to foams: the ratio is generally between 0.1 and 0.3 
for polymer melts and solutions [37l [38] . 

A measure of shape anisotropy for a given cell can be 
derived from the interface tensor [STJ [39] 
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where S refers to the surface of the cell and V its volume. 
The magnitude (second invariant) of g^j, a dimensionless 
scalar Q defined by 
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provides a simple gauge of shape anisotropy for an indi- 
vidual cell. Figure [l] shows the variation of each cell's 
Q in a foam at two different strains; Q typically ranges 
between and 1. Note that large Q values are typical in 
undeformed foams where the macroscopic stress is essen- 
tially zero, which indicates the individual cells are quite 
deformed. When the macroscopic stress is large, the cells 
are deformed and aligned. 

We consider the cell anisotropy changes during shear 
deformation of a disordered monodisperse foam with 512 
cells. The increment AQ for each cell is computed at 
elastic strain steps and across Tl cascades. Cells that 
become more isotropic (AQ < 0) experience stress re- 
lease. For elastic strain steps, the release domain^ the 
set of all cells with AQ < 0, typically contains between 
3 and 10 cells (Figure |4j left). The geometric transfor- 
mation at a Tl is more complicated. The set of cells 
that change topology across a Tl, the plastic zone^ can 



contain between 4 and 150 cells (Figure [4j center). The 
discrete changes in cell shape during Tl cascades signifi- 
cantly redistribute the cell level stress. The stress release 
domain across a Tl cascade contains considerably more 
cells (typically around 200), spanning the foam in all di- 
rections (Figure [4| right). The complex morphology of 
this set of cells is a percolating network-like structure. 

Network-like geometry of a surface can be quantified by 
the Euler characteristic x [40 ^ 41 , a topological invariant 
of a surface related to the number of holes and handles 
in the surface. It is equivalent to 2 — 2^, where g is the 
genus of the surface. For example, x = 2 for a sphere 
and X = for a torus. The invariant x of ^ surface S can 
be computed as 
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— / Kds, 
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where K is the Gaussian curvature. For structures result- 
ing from random spatial processes, large negative values 
of X are indicative of a network-like structure [42). 

Figure [5] shows the distribution of x of the stress re- 
lease domains for a shear simulation with 512 cells. The 
values of x form two distinct groups: those at an elastic 
strain step (purple) and those across a topological tran- 
sition (red). The x values at elastic strain steps follow 
a relatively smooth curve that peaks prior to the transi- 
tion to the stress plateau. In an isotropic foam, the prin- 
ciple eigenvectors of the cells are randomly distributed. 
When the foam is sheared, the cells begin to align along 
the principle eigenvector of the macroscopic stress. The 
observed peak in x signifies the progression of the cells 
towards alignment, and is a signature of large elastic de- 
formation and reorientation of the microstructure. 

Once the shear stress reaches a plateau (7 ^ 0.5), x 
is negative at, and only at, Tl cascades. The distinction 
between elastic deformations and Tl cascades is again 
highlighted in Figure [6] where each drop discontinuity 
in shear stress coincides with negative x and a Tl cas- 
cade. These correlations confirm the complex network- 
like structure of the stress release domains across Tl cas- 
cades, such as the example shown in Figure [T] The com- 
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FIG. 7. (Color online) (Left) The surface which bounds the 
network-like stress release domain of across a Tl cascade is 
shown in green, and the medial axis representation of this do- 
main is shown in purple. Here x — —70. (Right) The medial 
axis representation of the surface clearly reveals its network- 
like nature. The medial axis was computed by distance or- 
dering homotopic thinning of a voxelized image [42l [45] . 



FIG. 5. (Color online) The distribution of x of stress 
release domains for a foam with 512 cells sheared to a strain 
of 1.75. The black line shows the shear stress (^(7)), which 
reaches a plateau at 7 ~ 0.5. The purple and red points show 
the Euler characteristic x of the stress release domains at each 
elastic strain step and topological transition respectively, x 
forms two distinct groups: in the stress plateau, x > ^it 
elastic deformations, and x < at Tl cascades. 
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FIG. 6. (Color online) (Left) x versus the change in stress 
At at each elastic strain step (purple, both positive) and Tl 
cascade (red, both negative) for 7 > 0.5. (Right) Ar versus 
the size of the plastic zone approximated here by cells 
that change number of vertices, edges or faces, for 7 > 0.5. 
Drops in stress coincide with Tl cascades, and hence negative 
X indicates Tl cascades. 



bination of these results is a striking demonstration of the 
highly non-local effects of irreversible Tl cascades. The 
observation that a Tl cascade releases stress along a con- 
nected network spanning the foam is reminiscent of force 
chains and force networks in granular media [43\ 144] . 

Since simulations of significantly larger systems are not 
currently feasible with the Surface Evolver, the question 
of a long-term decay of the stress release domain cannot 
be addressed here. However, it appears likely that for 
cell numbers typical for experimental foam studies, the 
stress release reverberates through the entire system. 



In conclusion, we have demonstrated that Surface 
Evolver simulations of quasistatic shearing flow can pro- 
vide detailed insight into the processes that lead to 
elastic-plastic behavior of random soap froth. Our key 
observation is that the stress release associated with plas- 
tic flow events reverberates macroscopically throughout 
the foam sample, along a percolating network-like stress 
release domain. This highlights the complex rheology of 
foam, where discrete local events cause global geometric 
changes along with global redistribution of stress. 
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